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1 Introduction 



The approach to solid-body rotation of a fluid inside a rotating boundary is a familiar 
phenomenon with many applications. For instance, not only can we directly observe 
this phenomenon in the laboratory, but it may also play an important role in solar 
models, neutron stars and other environments. Greenspan & Howard (1963) give a 
fundamental analysis of the linearized version of this problem by considering a rotating 
axisymmetric container filled with a viscous incompressible fluid. They examine the 
behaviour of the fluid after the angular velocity of the container is suddenly changed 
by a small amount. Their solution consists of three distinct, time-separated phases: 
boundary layer formation, Ekman pumping and viscous relaxation. The bulk of the 
fluid spin-up (or down) occurs through Ekman pumping. Subsequent studies of the 
effect a stratification in density has on the Ekman pumping process (Walin 1969; 
Sakurai 1969; Buzyna & Veronis 1971; Hyun 1983; Spence, Foster & Davies 1992) 
found that some regions of the interior do not reach the final angular velocity until 
the viscous diffusion time has elapsed. The basic three stage structure of spin-up, 
however, remains intact where the intermediate Ekman pumping time results in a 
quasi-steady state in the fluid interior (see the review of Benton & Clark 1974 for a 
complete discussion of this subject). 

The spin-up of a compressible fluid has been studied in the context of rapidly 
rotating gases (Sakurai & Matsuda 1974; Bark, Meijer & Cohen 1978; Lindblad, 
Bark, & Zahrai 1994). These studies concluded that the basic structure for spin-up 
remained unchanged, but with some important and interesting differences. Numerical 
studies have also been carried out for both the linear and non-linear regimes (Hyun 
& Park 1990, 1992; Park & Hyun 1994) with the results lending support to the 
analytical results. The compressible fluid research, however, was motivated by a 
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desire to understand the dynamics of gas centrifuges where the effects of gravity are 
negligible compared to centrifugal forces. A density gradient results, but in a direction 
perpendicular to the angular velocity. 

In this investigation, we are primarily concerned with the possible astrophysical 
applications of the theory of stratified, compressible rotating fluids. Indeed one of the 
primary motivations for studying spin-up of a stratified fluid was to understand the 
sun's rotation (Howard, Moore & Spiegel 1967; Sakurai, Clark & Clark 1971; Clark, 
Clark, Thomas & Lee 1971). Though the sun does not have a solid outer crust, Ekman 
suction may arise because of solar wind torque, or boundary layer flow taking place 
between fluid interfaces. Ekman pumping in multilayer fluids was investigated by 
Pedlosky (1967) with his theory tested experimentally by Linden & van Heijst (1984) 
and O'Donnell & Linden (1992). Recently, the possibility that Ekman pumping 
may play a role in the synchronization of some binary stars has also been discussed 
(Tassoul & Tassoul 1990, 1992; Rieutord 1992). We, however, are motivated by 
the phenomenon of pulsar "glitches," a sudden slight increase in a neutron star's 
rotational frequency, and the resultant response of the fluid core. The simplified 
model of a rotating neutron star we consider includes the effect of compressibility 
of the inner fluid as well as a density gradient parallel to the gravitational field 
and angular velocity. A strong magnetic field, as exists on the surface, may alter 
the internal dynamics, but the presence of a magnetic field in the core is not well 
established and we choose to ignore it. The effects of a magnetic field which threads 
the core are discussed in greater detail in a forthcoming paper. 

2 Ekman pumping 

2.1 Fluid dynamics 



2 



To investigate the response of the fluid in a rotating container, we examine the 
usual, simple model of a cylinder of height 2L* and radius r c * rotating with angular 
velocity f2* (here and elsewhere an asterisk subscript indicates a dimensional variable 
or operator; quantities without this subscript are dimensionless.) When the angular 
velocity of the container is abruptly changed by a small amount, the differential 
rotation between the fluid and the top and bottom of the cylinder generates the 
"Ekman pumping" process. Unlike previous studies, we introduce an equation of 
state relating the mass-energy density to the pressure and to the composition. 
For a fluid with a kinematic viscosity v*, the Navier-Stokes equations of motion in a 
frame rotating with angular velocity 17* are 

P* + v, • V*v, + 2(0, x v*)^ = -V*p,+p*g* + ^p*V*fi^ + /o*i/*V2v*, (1) 

where g* is, 

8 * \ +g*e z z* < ' [Z) 
r\ is the cylindrical radius, and we take fi* = f2*e 2 with e z the unit vector in the 
z-direction. From this point onward we consider only the upper half-plane, noting 
that all quantities are symmetric about = 0. 

This particular form of g* represents a constant inward pointing acceleration, 
even though in a self-gravitating body the acceleration is a decreasing function of 
the radius, becoming zero at the center. However, as we will show below, for the 
parameter ranges applicable to a neutron star the significant dynamics occurs in a 
thin layer near the boundaries where gradients to the gravitational acceleration are 
negligible. 

As long as r c * is not too large the centrifugal acceleration is small compared to the 
gravitational acceleration and can be neglected. More precisely, we assume that finite 
Froude number effects can be ignored, i.e. F = 4£l1r c */ g* <^ 1, where F is the Froude 
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number. This results in a state of rotational equilibrium where the pressure p s * and 
density p s * are functions only of z*. The Navier-Stokes equation for the equilibrium 
system is 

9 (*\ 

We now look at a perturbed system in which the angular velocity of the boundary 
is suddenly changed by a relatively small amount Afi*. The resulting pressure and 
density are 

P* = Ps*{z*) + 5p*(r*,z*,U) (4) 
p* = Ps*{z*) + 5p*( r *> z*,t*)- (5) 

To first order in 5p*, ^p* and i>*(r*, z*, t*) we have 

<9v X X 

— + 2fi*e 2 x v* = V*5p* 5p* g* + i/*V*v*. (6) 

CC* p s * p s * 

We non-dimensionalize the equations by writing variables and operators as a dimen- 
sional constant times a non-dimensional variable or operator as follows, 

= (L*Afi*)v 

t, = (£ 1/2 2O0 _1 * 
= L* r 

Z;); == Z fig 

5p* = (2fi*p *-^*Afi*)5p 
^p* = (2Q,*po*L*AQ*/g*)5p 

Ps* = Po* Ps 
V, = (1/L,)V 
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where po* is a fiducial value for the equilibrium density. We also introduce the di- 
mensionless viscosity, or Ekman number, 



E 



(7) 



The Navier-Stokes equation for the perturbations is now 



E 1/2 — + e z x v = V5p 5pe z + £V 2 v, 




Ps Ps 



or in terms of the individual cylindrical components, 




(11) 



(10) 



(9) 



where (u, v, w) are the velocities in the (r, 6, z) directions. We need two more equa- 
tions in order to complete the formulation of the problem, an equation of state and 
the continuity equation. 

We describe the fluid in terms of the pressure and the concentrations of its con- 
stituent elements. Within the context of neutron stars these elements are mainly 
electrons, protons and neutrons. The equation of state then relates the density to 
these quantities, p* = p*(p*,Yi) where Yi is the concentration of the z-th particle 



The nature of the restoring force and the corresponding Brunt -Vaisala frequency 

is most readily calculated in the Lagrangian, as opposed to the Eulerian, formulation 

^dn the core of an equilibrium neutron star the Y$ are the concentrations that minimize the free 
energy through nuclear and weak interaction reactions. In the perturbations considered here, the 
fluctuation time scales are short compared to those for the weak interactions to adjust the ratio of 
neutrons to proton. The values of the Yi can thus be considered as fixed properties of the matter. 
If the equilibrium values of Yi give a stable stratification, buoyant forces will cause perturbations to 
oscillate with the Brunt -Vaisala frequency. 
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of the perturbations. We use <5g* for an Eulerian perturbation of a quantity q*, the 
difference between the actual and non-perturbed values of that quantity at a given 
point in space and time. A Lagrangian perturbation Aq* describes the change from 
the non-perturbed value an element of fluid experiences as it travels from one point 
to another. The two perturbations are related by a displacement vector field, 

Ag* = 5q* + ■ V*g *, (12) 

where (?o*( r ) is the non-perturbed quantity. The displacement vector field is related 
to the velocity by 

v, = (13) 
In non-dimensional notation ^ = (L#AQ*/2Q*) £ and 

v = E^y. (14) 

To relate the density and pressure perturbations, consider a fluid displacement in 
which some quantity Y is held constant i.e. AY = 0. The Lagrangian perturbations 
Ap* and Ap* are then related by 



Ap* = ( ^) Ap* = J-Ap„. (15) 

v OP* ) Y C Y* 



If the fluid is displaced adiabatically so that the entropy and composition are fixed, 
then cy* is the usual sound speed. We characterize the equilibrium relationship 
between the density and the pressure by 

dp s */dz* ( dp A _ 1 



dp s Jdz* \dp*J cq c^/ 
With fll2l)-( ll6|) we can relate <5p* and 8p*: 



(16) 



5p* = Ap* - £^ ■ Vp s * (17) 
= -^-5p* + f-j j p s *g*iz*, (18) 



with £j„ the ^-component of Once again, non-dimensionalizing we obtain, 

gp = ^J±\ 5p + fN?\ ^ _ Ky5p + N2p ^ (ig) 

where the Brunt- Vaisala frequency N* is 

N^d(^-P), (20) 

\ c eq* C Y*/ 

and the two dimensionless parameters Ky = g^L^/cy^ and N = N if /2Q :¥ are the 
"constant-V compressibility" and the normalized Brunt- Vaisala frequency, respec- 
tively In previous studies Ky was assumed to be negligible, but in self gravitating 
astronomical bodies Ky can be of order unity or much larger. Returning to our ex- 
ample of the neutron star, for instance, we can estimate the size of Ky. Using the 
values g* ~ 10 14 cm/sec 2 , « 10 6 cm and cy* ~ 10 9 cm/sec (Epstein 1988), we obtain 
Ky ~ 10 2 . N characterizes the influence of density stratification on Ekman pumping. 

At this point, it is worth noting that fll8|) is a generalization of more familiar 
formulas for low frequency density perturbations in a compressible stratified fluid. 
In particular, in studies of the terrestrial atmosphere one usually assumes that the 
fluid motions are slow enough to allow the pressure of displaced fluid to adjust to its 
surroundings so that 5p = and the first term in the left hand side of (jl8|) vanishes. 
This is equivalent to what we find in §6.4 of Pedlosky (1979). Furthermore, if we 
consider the incompressible limit, c Yjf — ► oo, then 

Nl - (21) 

which is the familiar incompressible formula for the Brunt-Vaisala frequency. 
The final equation is the continuity equation for the perturbations, 

^ = -V,.*. (22) 

Ps* 

— V r * • £ r * . (23) 
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With (|12D, the continuity equation becomes, 

5p* + £z*-ft- + Ps*V r * ■ £ r * + Ps*^ 1 = 0. (24) 
oz* oz* 

Using ( |T6| ) and (|I~8D, and taking the time derivative of (0), we get 

dp* =- w * + v r* ■ u * + — = 0. (25) 



Ps* c Y* ^* C Y* @ Z 

In non-dimensionalized units this is 



E 1 / 2 n 2 — ^-K Y w + — + - — (ru) = 0, 26 
at p s oz r or 

where we have introduced the dimensionless angular velocity, 

n = (27) 

We can now rearrange the complete set of perturbation equations in a more con- 
venient form. We use (plj) to eliminate Sp in dTT|) and take the time derivative to 
obtain, 



at 2 at p s at p s at 

where the "equilibrium compressibility" is 



* eq = (29) 

^eq* 

Since the constant-!" and equilibrium compressibilities are comparable, we write 

N 2 L L 

AK = K cq -K Y = = FN 2 — <1. (30) 

y* ''c* 

The variable 5p only occurs in the combination 5p/p s , so we define SP = 5p/p s - The 
final equations are now 

£l/ i«-» = -^ +£ ( v2 -^)" < 31 > 



S 



E y^v + u = Eh 2 -\)v (32) 



dt " V " r 2 

E^ = -N 2 w - E^~SP + E^Ak^-SP + £ 3/2 V 2 |^. (33) 
ot 2 oz at at at 

E 1 / 2 n 2 ^-5P-KYW + -^-w + --^-(ru) = 0. (34) 
at oz r or 



Note the harmonic restoring force provided by the Brunt- Vaisala term in (|33|). 
The above four equations, (]3T|)-(|34|), describe the evolution of the four unknowns, 
v and SP. We have dimensionless parameters, E, N 2 , An, Ky and Q 2 . In order to 
reduce the parameter space we consider only slow rotation, fi 2 < 1, and Ak <C 1. 
Both n 2 and A k are easily included in the general solution, but since they only 
appear as a product with E 1 / 2 their effects are small and we do not consider these 
terms in what follows. Table 1 contains definitions of the dimensionless parameters 
and representative values for a neutron star. 

The presence of the compressibility Ky distinguishes this set of equations from 
earlier studies (Walin 1969; Sakurai 1969; Clark et al. 1971). Previous studies chose 
to emphasize the effects of temperature on the density of the fluid. Specifically, the 
density was considered a function of the temperature and the stratification was a 
result of a temperature gradient which was imposed by the boundary conditions. 
The dynamical significance of the stratification and Brunt-Vaisala frequency arose 
through the effects of temperature diffusion and the heat equation. This approach 
is not appropriate to the astrophysical cases with which we are primarily concerned. 
In neutron stars, for example, thermal effects have a negligible result on the fluid 
dynamics, whereas compressibility is quite significant. We, therefore, focus on the 
dependence on the equation of state. 

2.2 Boundary values and initial conditions 
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To obtain a unique solution to (|31])-(|34|), we need to specify both the boundary 
and initial conditions to our problem. There are, in essence, two approaches to 
take at this point. The most complete method is to state that initially the fluid 
rotates uniformly with the cylinder, and solve for the behaviour of the fluid after the 
angular velocity of the cylinder changes with the Laplace transformation technique 
(Greenspan & Howard 1963). A simpler, and more physically elucidating approach, 
although less rigorous, used by other researchers in the field (Walin 1969; Sakurai 
1969; Barcilon & Pedlosky 1967) entails recognizing that different physical processes 
take place on widely different time scales in different regions of the fluid. We will 
follow this latter approach. 

If the Ekman number E, or dimensionless viscosity, is sufficiently small, the be- 
haviour of the fluid following an abrupt change in rotation rate of the container can be 
viewed as three distinct physical processes which occur on time scales fi" 1 , E~ l l 2 VL~ l 
and E^VL' 1 . The most rapid process is the formation of a viscous boundary layer. 
Following the impulsive change of rotation of the cylinder, a viscous Rayleigh shear 
layer forms on the upper and lower surfaces in a time scale on the order of a rotation 
time (tfo ~ ^^T 1 )- Within this region the gradient in the azimuthal velocity results 
in an imbalance between the centrifugal and pressure gradient forces causing fluid 
to flow radially. This radial flow in the boundary layer establishes a secondary flow 
where fluid in the interior is pulled into the boundary layer to replace the flow in the 
Ekman layer, creating an opposing radial flow in the interior fluid that satisfies conti- 
nuity requirements. This Ekman pumping spins the interior of the fluid up in a time 
scale of order .E -1 / 2 ^" 1 . With our choice of dimensionless variables this corresponds 
to a dimensionless time, £g ~ 1. Finally, residual oscillations decay in the viscous 
diffusion time t m ~ E^ 1 ^ 1 . 

Since the principal goal of this investigation is to understand the effects of the 
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stratification and compressibility on the Ekman pumping in the interior of the fluid, 
we expand (|3~T|)-(|3~i|) in powers of E 1 / 2 and isolate the equations relating to Ekman 
pumping. The initial velocity distribution for the Ekman pumping equation is equiv- 
alent to the final velocity distribution of the boundary layer which forms during the 
first phase. Following Walin (1969), we formulate the boundary condition in terms 
of the continuity of the velocity perpendicular to the Ekman boundary just outside 
the boundary layer, 

£1/2 

w(z = ±l) = T-^-(Vxv) 2 (35) 
E 1/2 1 d 

= T^r-!-M. (36) 

It is critical to the dynamics of Ekman pumping that the vertical velocity at the 
boundary layer is 0(E 1 ^ 2 ). This standard result (see, e.g. Pedlosky 1979) can be 
understood by scaling arguments. The imbalance between the centrifugal forces and 
pressure gradient forces in the boundary layer drives the Ekman pumping process. 
The thickness A of the boundary layer is 0(E 1 ^ 2 ) since the viscous terms in the 
dimensionless Navier-Stokes equation is EV 2 ~ E/X 2 = 0(1). The mass flux within 
the boundary layer is M\ oc A = 0(E 1 ^ 2 ). The net mass flux M z oc w perpendicular to 
the boundary layer is of the same order as M\ giving w = 0(E 1 ' 2 ). The sidewalls also 
have an 0(E 1 ^ 2 ) boundary layer. As shown in Pedlosky (1967) the vertical velocity 
through this layer is inhibited by buoyancy forces and is 0(E 1 ^ 2 ) leading to a net 
mass flux of 0(E). The boundary condition of the radial velocity is then u = 0(E) 
at r c . Since we keep terms only up to 0(E 1 ^ 2 ), as described in the following section, 
this gives u(r c ) = 0. 
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3 Solutions 



To solve fl3"T|)-(|3"4T) perturbatively we expand each fluid variable q as a series q = 

go + E 1 ^ 2 qi + Eq2 + ■ ■ ■■ Collecting terms of a given power of E lj/2 , we obtain a set of 
equations governing each order in the expansion. 
We find that the 0(1) equations are, 

v = ^-6P (37) 

u = (38) 

w = (39) 

did.. . . 

- K Y w + —w + -—(ru ) = 0, (40) 
Oz r Or 



and the 0(E 1 ^ 2 ) equations are, 

7,1 = 

Or 



vi = -i-SP, (41) 



= -ui (42) 

Kywi + —wi + -T-(r Ul ) = 0. (44) 
oz r or 



We define 4> = —d5P /dt, so that (|37|) and (03) become 



<9r 

and Q4"3] ) and ([P]) are now, respectively, 



wi = ^-0 (45) 



iVV = (46) 
oz 

-K Y ( d \ d 1 ( d \ 19.9, n 

^ + ^ IteV + rd?^ = °- (4?) 
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Assuming N 2 varies slowly over z, we treat it as a constant and simplify ( |47D to 

az z (72 r or \ or J 

After taking the time derivative, the boundary condition, (1351) , is 



d d , iV 2 l5 / 9 \ 

±-=— — r— , at* = ±l. (49) 



By setting = Z(z) i?(r) T(t), flilf ) becomes, 



1 d 2 _ Kyrf iV 2 1 d ( d \ n 
Zdz^ Z -YTz Z+ ll-rTr{ r Tr R )= - (50) 

The solutions to the spatial functions are 

Z = AeP+ z + Be 15 - 2 (51) 
R=J (kr), (52) 

where 

/?± = i(Ky±(Ky 2 + 4fc 2 JV 2 )*) (53) 
The symmetry of the boundary condition, W\{z = 0) = 0, relates A and B: 

A = -P—B. (54) 
P+ 

The constant B is arbitrary, and we choose it so that Z(l) = 1. This leads to 

B = (e?- - ^ef^J 1 . (55) 
The possible values of k are determined by the boundary condition at the sidewall at 



r = r c , i.e. Ui(r c ) = 0. From fl45|) and (p2|) we see that this condition corresponds to 
Ji(k m r c ) = for m = 0, 1, 2, . . .. The first zeros of J\ are fc m r c = 0, 3.8317, 7.0156, . . .. 
The solution k m = has v = everywhere and is of no interest. 
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We utilize the boundary condition to determine the time dependence of 0. Putting 
our solution for R and Z into (0), we obtain the differential equation, 

-k 2 N 2 Z(l) . , 

T T = 7= H7,\ T, 56 

dt ^2 f (1) 1 1 

whose solution is, 

-Ult 



T(t) = e _wt , (57) 



where, 



h 2 N 2 -1 

uj = ^-(Af3 + e^ + B(3^-y , (58) 

with A and B defined as above. 

There are two interesting limiting cases. The first is that of no stratification 
iV — > 0; the second is that of an incompressible fluid Ky — > 0. Let us consider the 
first of these which gives, 

/ k 2 N 2 \ 

*-"*( 1 + -^J (59) 

k 2 N 2 

(3. « (60) 
A; 2 iV 2 

S«l + ^2_(l + e ^/« y ) (61) 

Ky 

k 2 N 2 

A « (62) 

and, 

w « -=-^ . (63) 

v/2(e^-l) V ; 

Equation d63f) shows that for large compressibilities the Ekman spin-up time scale u^ 1 
grows exponentially with Ky. The second limiting case, Ky — > 0, gives (3± ~ ±k m N 
and 

k m Ncoshk m N 

uj Pd — p . (64) 

V2 sinh fc m iV 

This matches the Ky = 0, iV ^ solution which was obtained by Walin (1969). 
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We are now in a position to write the complete solution for the quantities and 

oo 

= ]T C m Z m {z) Jo(fc m r)e-^*. (65) 

m=l 

The velocity, t>o, is found from the relationship, 

y = -^vo, (66) 

which gives, 

oo u 

v (r,z,t) = - £ —CmZ^Mk^e-^ + v^z). (67) 

m=l ^ 

The last term represents the final velocity due to Ekman pumping. If we take the 
frame of reference as that rotating with the cylinder before impulsive spin-up, the 
final velocity at the boundary of the interior fluid is 

v^ir.z = ±1) = r. (68) 

We determine C rn from the initial state of the fluid, 

v (r, z = 1, t = 0) = = - £ -^^Jiikmr) + r. (69) 

m=l U;m 

The coefficients, given by the standard equation for a Fourier-Bessel series, are 

h C 9 

nj m\ J m ^ (70) 

UJ m k m J2\km) 

and the final velocity is 



Vooir.z) = 2 ^ i T 77 , Ji(k m r)Z m (z), (71) 

m= l K m J2{K m ) 



where we have chosen r c — 1 (r c * = L*). 
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4 Discussion 



The time dependence of the Ekman pumping process is exponential with characteristic 
time 1/uj. We plot the value of u as a function of k m N in figure 1 for different values 
of the parameter Ky. Larger N, corresponding to greater density stratification gives 
larger uj and reduced characteristic time. That is, a strongly stratified fluid spins up 
much quicker than a non-stratified fluid. On the other hand, an increased value of 
the compressibility Ky slows the pumping process for a given k m N . The spin-up time 
uj^ 1 decreases with increased stratification because stratification isolates much of the 
fluid from the pumping process. 

Figures 2 and 3 show that the rotation state at the end of the Ekman pumping 
stage is not that of a solid body. The ordinate Z(z) is proportional to the final 
azimuthal velocity, with Z = 1 being the largest possible spin-up. Larger values of 
k m N leave more of the internal fluid unaffected by the Ekman pumping process. In 
contrast, in a homogeneous fluid, N — 0, Ekman pumping brings the entire fluid to 
an angular velocity equal to that of the boundary. The compressibility Ky further 
decreases the amount of pumped fluid, as we can see by comparing figure 2, for 
Ky = 0, with figure 3, for Ky = 10. 

Compressibility thus decreases the efficacy of Ekman pumping both by lengthening 
the spin-up time and by decreasing the amount of affected fluid. To convey a clearer 
picture of how strong the effect of Ky is, we plot in figure 4 the final angular velocity 
of the fluid at its central (z = 0) layer as a function of k m N for different values of Ky. 
Though there is little change between Ky — and Ky — 1, the internal final angular 
velocity is strongly suppressed as Ky increases to 10. 

We point out that for canonical values N 2 6 and Ky ~ 10 2 for a neutron star, 
the thickness of the layer affected by Ekman pumping is much smaller than the radius 
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of the cylinder. Gradients in the gravitational acceleration are therefore small in these 
layers, justifying our original assumption of constant g* in (2). 

In figure 5 we plot the average spin- up of the fluid (Z) as a function of the 
normalized Brunt-Vaisala frequency N for the two lowest order modes, k\ and A)2- 
We see that even modest values of N prevent most of the fluid from spinning up during 
the Ekman pumping phase. The state of the fluid after a time scale of « E^ 1 / 2 ^ 1 
is, thus, one of non-uniform rotation. The process of viscous diffusion, which operates 
in a time t m ~ E^ 1 ^ 1 eventually brings the fluid into solid-body rotation. 

The case of spherical geometry was studied by Clark et al. (1971), where they 
found that the solution for a sphere is qualitatively similar to that of the cylinder. 
That is, the final state of non-uniform rotation also exists in the sphere, but the 
geometry of the layer that gets Ekman pumped is modified. 

We find a particularly interesting application of these phenomena is the response 
of the interior of a rotating neutron star to a glitch, a sudden small change in the 
rotational velocity. Within the star there exists a significant stratification due to the 
strong gravitational field and the equilibrium concentrations of protons, neutrons and 
electrons. Reisenegger & Goldreich (1992) estimated a value of A* « 500s _1 for a 
neutron star. For a canonical value of fi* 100s" 1 , we obtain N as 2.5. This is large 
enough to have a significant effect on the length of time the core of the star needs to 
come into rotational equilibrium. We explore these issues in a forthcoming paper. 

We would like to thank Angela Olinto for her valuable input and advice. M.A. 
was supported in part by the DOE at Chicago, by NASA grant NAG 5-2788 at Fermi 
National Laboratory, and through a collaborative research grant from IGPP/LANL. 
This work was carried out under the auspices of the U.S. Department of Energy. 
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Figure Captions 

Figure 1. The spin-up characteristic time, w, as a function of k m N for varying 
values of n Y . 

Figure 2. The final azimuthal velocity as a function of depth for an arbitrary 
value of the radius. A value of Z — 1 is complete spin-up, while Z = is no spin-up, 
with ky = 0. 

Figure 3. As in figure 2, but with k y = 10. 

Figure 4. The final velocity of the central layer of the fluid [z = 0), as a function 
of k m N, for different values of k y . 

Figure 5. The average final spin-up of the fluid as a function of the stratification. 
The top two curves (solid and dashed lines) were calculated for an incompressible fluid, 
ky = 0, while for the bottom two curves (dotted and dash-dotted lines) ky = 10. 
With a highly compressible fluid (/ty = 10) even very small values of iV result in very 
little spin-up from Ekman pumping. 
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parameter 


formula 


value 


E 


vJ(2n*Ll) 

T 1 \ -T- J 


10- 7 


N 2 




6 


K Y 




10 2 
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^lr c */g* 


10" 4 


An 


N*L*/g* 


10" 4 






1CT 2 



Table 1. The dimensionless parameters. The values quoted are order of magni- 
tude estimates for a characteristic pulsar. 
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